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Abstract. The eigenfunctions e l ^ ,x ' of the Laplacian on a flat torus have uniformly 
bounded LP norms. In this article, we prove that for every other quantum integrable 
Laplacian, the LP norms of the joint eigenfuntions blow up at least at the rate lly'fcllip > 



C(e)A fc 4p when p > 2. This gives a quantitative refinement of our recent result }TZ1| 
that some sequence of eigenfunctions must blow up in LP unless (M,g) is flat. The better 
result in this paper is based on mass estimates of eigenfunctions near singular leaves of the 
Liouville foliation. 



0. Introduction 



This paper, a companion to ||TZ1|| , is concerned with the growth rate of the LP- norms of 
L 2 -normalized A-eigenf unctions 



A(fj = XjCfj, ((fj, (p k ) = 5 jk 

on compact Riemannian manifolds (M, g) with completely integrable geodesic flow G t on 
S*M. The motivating problem is to relate sizes of eigenfunctions to dynamical properties 
of its geodesic flow G l on S*M. In general this is an intractable problem, but much can be 
understood by studying it in the framework of integrable systems. To be precise, we assume 
that A is quantum completely integrable or QCI in the sense that there exist P 1; . . . , P n G 
\I/ 1 (M) (n = dim M) satisfying [Pi, Pj] = and such that their symbols (jq, . . . ,p n ) satisfy 
the standard assumption that dpi A dp2 A ■ ■ ■ A dp n ^ on a dense open set Q C T*M — 
whose complement is contained in a hypersurface. Since {pi,Pj} = 0, the pi, . . . ,p n generate 
a Hamiltonian M n -action given by the joint Hamiltonian flow £) := exp(tiX pi ) o ... o 
exp(tiX Pn )(x, £), where for j = 1, ...,n, X p . denotes the Hamiltonian vector field of pj. The 
associated moment map V : T*M — > M. n is V = (pi, ■■-,p n ) and the orbits of this action foliate 
T*M - and since S*M is preserved they also foliate S*M. See [fTZl| , |TZ2| for a list of 



well-known examples. We will refer to this foliation as the Liouville foliation. Throughout, 
we will assume that the leaves of this foliation satisfy Eliasson's non-degeneracy hypotheses 
(see 



[FJJ and Definition [T|) . 
The main result of [[TZ1|| was that the L°°-norms of the L 2 -normalized joint eigenfunctions 
{(fix} of (Px,...,P n ) are unbounded unless (M,g) is flat. In this paper, we investigate 
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the blow-up rates of ZAnorms of sequences ||<£>a||l*> under the Eliasson non- degeneracy 
assumption. 

Theorem 1. Suppose that (M,g) is a compact Riemannian manifold with completely inte- 
grate geodesic flow satisfying Eliasson's non- degeneracy condition. Then, unless (M,g) is a 
flat torus, there exists for every e > 0, a sequence of eigenf unctions satisfying: 

<Pk\\L- > C{e)\\' e . 

Vk\\u> > C{e)xf~\ 2 <p< oo. 

This result is sharp in the setting of all completely integrable systems. It is based on the 
existence of a codimension one singular leaf of the Liouville foliation. Here, a singular leaf 
is an orbit of the R™ action with dimension < n, see Definition (lip. When codimension I 

leaves occur, the estimate becomes ||v?£:||l°° > C(e)\ k 
Existence of a singular leaf is given in the following: 

Lemma 2. Suppose that (M, g) is a compact Riemannian manifold with completely integrable 
geodesic flow. Then the Liouville foliation of (M, g) contains a singular leaf unless (M, g) is 
a flat torus. 

The proof only shows that some singular leaf occurs, but does not determine its codi- 
mension. As will be discussed in §[[[ it is quite plausible that codimension n — 1 singular 
leaves often occur for Hamiltonian M. n actions on cotangent bundles. Such leaves correspond 
to closed geodesies which are invariant under the M. n action. In such cases, the estimate 

improves to ||^fc||i,°o > C(e)X k 4 . This estimate agrees with a recent result of Donnelly 
|D|] on blow-up rates of L°° norms of eigenf unctions on compact (M, g) with isometric S 1 

n-l 

actions. He finds sequences of joint eigenfunctions which blow up at the rate A fc 4 ; they are 
precisely the same kind of eigenfunctions which are associated to singular leaves, as we now 
describe. 

The analytic ingredient of our proof gives an estimate on the growth rate of modes and 
quasimodes associated to a singular leaf of the Liouville foliation. It is here we assume that 
the singular leaf is Eliasson non-degenerate. 

Lemma 3. Let V~ 1 {c) with c = {c\, C2, c n ) e W 1 be a singular level of the Lagrangian 
fibration associated with the quantum completely integrable system Pi = a/A, P n . Suppose 
7 ?_1 (c) contains a single compact, Eliasson non- degenerate orbit A := W 1 ■ (vq) of dimension 
£ < n and that each connected component of it (A) C M is an embedded l-dimensional sub- 
manifold. Then, there exists a sequence tpk of L 2 -normalized joint eigenfunctions of P\, ...,P n 
with Pjipk = Ejkfk, Ejk = AfcCj + 0(1); j — 1, ...,n such that for any e > 0, 

n — £ 

We also have that 

(n-£)(p-2) c 

(ii) \\fk\\Lp > C(e)X k 4p , 2 < p < oo. 
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For simplicity we have stated Lemma (|3|) in the case where the level V~ l (c) contains a 
single Eliasson non-degenerate orbit, but the result extends in a straightforward fashion to 
the case of multiple singular orbits. 

Some final words to close the introduction. First, in |}rZ2|| we give a different approach 



to eigenfunction blow-up in the completely integrable case which in some ways is closer 
in spirit to the approach of this paper. Namely, we relate norms of modes to norms of 
"quasi-modes," i.e. approximate eigenf unctions associated to Bohr-Sommerfeld leaves of the 
Liouville foliation. In that paper, we also give a number of detailed examples of quantum 
completely integrable systems such as Liouville tori, surfaces of revolution, ellipsoids, tops 
and so on. 

Second, we would like to emphasize that results on completely integrable systems are rele- 
vant to quite general classes of Riemannian manifolds. This is because one can approximate 
any geodesic flow by a completely integrable flow (its Birkhoff normal form) near a closed 
orbit or invariant torus. On the quantum level, one may approximate the Laplacian by its 
quantum Birkhoff normal form. The quasi-mode analysis of this paper is therefore of some 
relevance to general Riemannian manifolds. What is special about quantum integrable sys- 
tems is the particularly close relations between modes and quasi-modes, which allows one 
to draw strong conclusions about the blow-up of sequences of modes from that of the corre- 
sponding quasi-modes. It is quite reasonable to suppose that the same kind of phenomenon 
occurs much more generally. In particular, we would conjecture that any compact (M, g) 
with a stable elliptic orbit has a sequence of eigenfunctions whose L°° norms blow up at the 

rate A fc 4 . There surely exist quasi-modes with this property, and the difficulty is to show 
that this implies the existence of modes with this blow-up rate. We hope to pursue this 
direction in the future. 

1. Geometry of completely integrable systems 

This section is devoted to the geometric aspects of our problem. We first prove Lemma 
(0) and then give background on Eliasson non- degeneracy of singular orbits. We begin with 
some preliminary background on completely integrable systems, in part repeated from |[LZ1 



As mentioned in the introduction, the moment map of a completely integrable system is 
defind by 

(1) V = ( Pl ,...,p n ):T*M^BcR n . 
The Hamiltonians pj generate the IR n -actin 

$ t = exp txE Pl o exp t 2 E P2 ■ ■ ■ o exp t n E Pn . 

We denote <3> r orbits by R n ■ (x, £). 

By the Liouville- Arnold theorem, the orbits of the joint flow $ t are diffeomorphic to 
M fe x T m for some (k,m), k + m < n. By the properness assumption on V, each connected 
component of a regular level is a Lagrangean torus, i.e. 

(2) V-\b) = AM (6) U • • ■ U A^\b), (b e B reg ) 
where each A^(fe) ~ T n is a Lagrangian torus. 
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We are particularly interested in singular orbits. To clarify the notation and terminology, 



we repeat the definition from ||TZ1|| : 

Definition: We say that: 

• b G B S i n g if is a singular level of the moment map, i.e. if there exists a point 
(x, £) G V^ib) with dpi A ■ • • A dp n (x,£,) = 0. Such a point (x, £) is called a singular 
point of V . 

• a connected component of "P _1 (6) (6 G B S i ng )is a singular component if it contains a 
singular point ; 

• an orbit ~R n ■ of <& t is singular if it is non-Lagrangean, i.e. has dimension < n; 

• b G -B re9 and that V~ l {b) is a regular level if all points G 7- >_1 (6) are regular, i.e. 
if dpi A • • • A dp n (x,C) ^ 0. 

• a component of ( 6 G £> S m 5 U -B re9 ) is regular if it contains no singular points. 

When b G B sing we first decompose 
(3) V-\b) = Vf 3=1 T^ ng {b) 

the singular level into connected components r*^L(&) and then decompose 



(4) r& ,(6) = UjUR n ■(**,£ 



A- J 



each component into orbits. Both decompositions can take a variety of forms. The regular 
components r*^ (6) must be Lagrangean tori by the properness assumption. A singular 
component consists of finitely many orbits by the Eliasson non-degeneracy assumption (see 
below). The orbit W 1 ■ (x, £) of a singular point is necessarily singular, hence has the form 
R x T m for some (k, m) with + m < n. Regular points may also occur on a singular 
component, whose orbits are Lagrangean and can take any one of the forms M fc x T m for 
some (k, m) with k + m = n. 

1.1. Proof of Lemma In this section, we prove Lemma (0): the Liouville foliation 
of metrics with completely integrable geodesic flows always contain singular leaves unless 
(M, g) is a flat torus. 

Proof. The hypothesis that the Liouville foliation is non-singular has two immediate geo- 
metric consequences: 

• (i) By Mane's theorem [|M| , (M,g) is a manifold without conjugate points; 

• (ii) The (homogeneous) moment map V : T*M — —>■ W 1 — is a torus fibration by T n . 

Statement (ii) follows from the Liouville- Arnold theorem. On a non-singular leaf, we must 
have dpi A • • • A dp„ ^ 0. Since this holds everywhere, V is a submersion; and since it is 
proper, it is a fibration. The fiber must be T n , again by the Liouville- Arnold theorem. Since 
V is homogeneous, the image V(T*M — 0) = 1R + ■ V(S*M). Since V is a submersion, the 
image is a smooth submanifold of 5 1 ™ -1 hence must be all of S^ 1 . 

By (i) it follows that M = M /V where (M, g) is the universal Riemannian cover of (M, g) 
(diffeomorphic to M n ), and where V = 7Ti(M) is the group of covering transformations. 
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We claim that (ii) implies tti(M) = Z n . Indeed, T*M — is a double fibration (with 
indicated fibers): 

T*M-0 

7T / (R n -0) (T n ) \ V 

M R™ - 

By the homotopy sequence of a fibration tt : E — > B, 

■ ■■TT q (F) -> ir q (E) -> 7r g (B) -> 7r ff _!(F) ► tt (£) -> tt (B) -> 

and using that 7r 2 (T n ) = 1 and that 7r 2 (M) = 1 by (i), we obtain 
(5) 

1 -> TTi^" 1 ) -> 7n(5*M) -> 7Ti(M) -> ^(.S"- 1 ) -> 7T (£*M) -> 7T (M) -> 1 

1 -> 7T 2 (5*M) -> 7T 2 (5 rl - 1 ) -> 7n(T") -> 7Ti(5*M) -> 7Tl (5' n— -> 7T (T n ) -> 7T (5*Af) -> ^(S^ 1 

Since 7r 2 (A^) = 7r 2 (A^) (where N) is the universal cover, and since S*M = S™-^ 1 x IR n , we 
have 7T2(S*M) = 7r 2 (5' n ~ 1 ). From its definition, we see that the homomorphism 7r 2 (S'*M) — > 
7r 2 (5' n_1 ) is an isomorphism, hence the second sequence simplifies to 

(6) 1 -> 7n(T") -> 7ri(5*M) -> miS"- 1 ) -> vr (T n ) -> tt (5*M) -> 7T (5 ,n - 1 -> 1 

Let us first assume that n > 3. Then tt\{S*M) = Ti\{M) and 7Ti(S' n_1 ) = 1, so we get 

1 — > ^l(T n ) — > 7Ti(M) — > 1, 

i.e. 7Ti (T n ) = tti(M). 

We now consider dimension 2. Actually, it follows by a classic result of Kozlov jK] that 
the only surfaces that can possibly have a completely integrable geodesic flow (even with 
singularities) are M = S 2 ,T 2 . We can easily disqualify M = S 2 since S*S 2 = IRP 3 does not 
fiber over S 1 . 

It follows then that M is diffeomorphic to R^/Z", i.e. M is a torus. Since g has no 
conjugate points, the proof is concluded by Burago-Ivanov's theorem pi . 



□ 

1.1.1. Remarks on the geometry. We would like to add two remarks on the geometry. We 
will not use them in the proof of our results, and therefore only sketch the issues. Our second 
remark could lead to a substantial improvement of our main result. 
Remark (i): On flat manifolds. 



In our previous article |[TZ1|| , we showed that eigenfunction blow-up occurs in the integrable 
setting unless (M, g) is flat. Here, we are showing that it occurs unless (M, g) is a flat torus. 
We briefly explain why Liouville foliations of other flat manifolds must have singular leaves. 

The Riemannian connection of (M, g) gives a distribution of horizontal planes in the unit 
sphere bundle, i.e. T(S*M) — H © V where V is the tangent bundle to the fibers of 
7T : S*M — > M. If g is a flat metric, then H is involutive and hence S*M is foliated by 
smooth n-manifolds L (which we will call horizontal leaves). (M,g) is also the Riemannian 
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quotient of a torus T n by a finite group G of isometries, and the derivative of the projection 
p : T n — > M takes the horizontal distribution of T n to that of M. Since the horizontal leaves 
of T n are compact Lagrangean tori, it follows that the horizontal leaves of (M, g) are compact 
Lagrangean submanifolds of T*M — 0. Thus, the geodesic flow of (M, g) leaves invariant a 
non-singular Lagrangean foliation. 

However, the horizontal leaves are not orbits of a Hamiltonian W n action in general. (We 
thank W. Goldman for several helpful discussions on this point). In fact, the leaves have the 
form L(G)\R"/Z ri where L(G) is the linear holonomy group (a finite group). So the leaves 
are quotients of tori, but L(G) is generally not a subgroup of T" so the leaves are finite 
quotients of tori. If L{G) acts freely on T*M — 0, then the leaves are tori, but there is no 
]R n action generating them. 

The natural 'Hamiltonians' on T*M — belong to the algebra I G where X =< I±, . . . ,I n > 
is the algebra generated by the action variables Ij(x,£) = £j on the covering torus, and where 
X G denotes the G-invariant elements. Since we require that the R n action be homogeneous, 
we only consider functions qj{h, ■ ■ ■ ,I n ) which are homogeneous of degree one relative to 
the natural M + action on T*M — 0. A completely integrable system on T*M — is given 
by a choice of n generators of this algebra. The M" action they generate will leave invariant 
the horizontal foliation, but the leaves of that foliation are not in general the orbits of the 
action. What the proof above shows is that any such action must have some singular orbits, 
of lower dimension than the leaf it is contained in. 

Remark (ii): On existence of high codimension singular leaves. 

The simple proof above on existence of singular leaves does not give any information 
on the codimension of such leaves. As mentioned in the introduction, our lower bounds 
llVftlU 00 — C(e)^fc improve as this codimension q increases. Can a better argument prove 
the existence of codimension n — 1 leaves in generic cases? 

Here is heuristic existence argument: the image C of T*M — under the moment map V 
of a homogeneous IR n -action is a convex polyhedral cone in IR n . In Lemma @, we observed 
that V cannot be a torus bundle over R n — unless (M, g) is a flat torus. In all other cases, 
C will have a boundary over which the singular orbits 'fiber'. 

The geometric problem is to determine the geometry of C. In the case of Hamiltonian 
torus actions, the moment polyhedra have been studied by Lerman et al. |I], L2| . They do 



not seem to have been studied for general Hamiltonian IR n actions. However, it is natural to 
conjecture that as long as C has a non-empty boundary, there will boundary faces of each 
codimension and in particular there will be one-dimensional edges. If this is true, then our 

lower bound improves to >> X k 4 

1.2. Eliasson non-degeneracy. Let c := (ci,...,c n ) € B S i ng be a singular value of the 
moment map, V. Suppose that: 

(7) rank (dpi, dp n ) = rank (dpi, dpk) = k < n 

at some point u G V^ x (c). Denote the orbit through v by 1R™ • (v ) := {exptiE pi o ■ ■ ■ o 
exp t fc S Pfc (f ); t = (ti, ...,t k ) G lR n }, which we henceforth assume is compact. Our aim is to 
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define the notion of Eliasson non-degeneracy of such a singular orbit; this is a straightforward 



extension of the notion defined in |E| for fixed points of symplectic maps. 

By the Liouville-Arnold theorem, the orbits of the joint flow $t are diffeomorphic to 
K>fc x T m for some (k,m), k + m < n. Let us first consider regular levels where m = n. By 
the properness assumption on V, a regular level has the form 

(8) V-\b) = A«(6) U • • • U A^\b), (b G B reg ) 

where each A® (6) ~ T n is an n-dimensional Lagrangian torus. The classical (or geometric) 
multiplicity function m c i(b) = i.e. the number of orbits on the level set is 

constant on connected components of B reg and the moment map V is a fibration over each 
component with fiber (§). In sufficiently small neighbourhoods Q( l >(b) of each component 
torus, A^(6), the Liouville-Arnold theorem also gives the existence of local action-angle 
variables (I®, In, &[ , ■ ■■■,Q < n) in terms of which the joint flow of S Pl , E Pn is linearized 
[AM]. For convenience, we henceforth normalize the action variables lf \ I n so that 7, — 



j 



0; j — 1, n on the torus Aw (6). 

Now let us consider singular levels of rank k. We first observe that dpi, ...,dpk (in the 
notation of (|7|)) are linearly independent everywhere on M n -(t>o). Indeed, the Liouville-Arnold 
theorem in the singular case [AM] states that there exists local canonical transformation 

ip = z/>(7, 9, x, y) : M 2n -> T* M - 0, 

where 

7 = (7i, ...,I k ),9 = (Ox, 9 k ) eR k , x = (xi, x n . k ),y = (y 1} ...,y n - k ) e M n - fc m 
defined in an invariant neighbourhood of M n • (f o) such that 

(9) pi o = Ii (i = l, k), 
and such that the symplectic form uj on T*M takes the form 

k n—k 

(10) iI)*uj = S2 dI J A dd i + A dy i- 

3=1 3=1 

As a consequence of the normal form (PH|), it follows that there exist constants with 
z = A; + 1, n and j = 1, k, such that at each point of the orbit, M. n ■ (vq), 

k 

(11) dpj = y^ j c ij dp j . 

3=1 

Since dpi, ...,dpk are linearly independent in a sufficiently neighbourhood U of v G 7- ,_1 (c), 
the action of the flows corresponding to the Hamilton vector fields, H Pl , H Pfc generates a 
symplectic M fc action on 7 :>_1 (co) H U. An application of Marsden-Weinstein reduction yields 
an open 2(n — /c)-dimensional symplectic manifold: 

(12) := V-\co) n U/R k , 

with the induced symplectic form, a. We will denote the canonical projection map by: 



7r k :V- L (c )nU 



k- 
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Since {pi,Pj} = for all i,j = 1, ...,n, it follows that Pk+i, ■ ■-,p n induce C°° functions on E*,, 
which we will, with some abuse of notation, continue to write as Pk+i, ■■■,p n - From (|lT|), it 
follows that 

dpi{7Tk{v )) = 0; i = k + 1, ...,n. 

Here, we denote the single point ^(R" • (v )) by 7ik(vo). 

To describe the Eliasson construction, let (M 2n , a) be any symplectic manifold and C, the 
Lie algebra of Poisson commuting functions, fi G C°°(M), with the property that: 

fi(m) = dfi(m) = 0; i = 1, ...n, 

for a fixed m G M. Let <2(2n) denote the Lie algebra of quadratic forms and consider the 
Lie algebra homomorphism, 7i : C — > Q(2n) given by 

H(f) :=d 2 f(m). 

The result of Eliasson says that, if TC(fi, / n ) is a Cartan subalgebra of <2(2n), there 
exists a locally-defined canonical mapping, «;:[/—> [To, from a neighbourhood, f/, ofmeM 
to a neighbourhood, f/ , of G M n , with the property that: 

(13) Vi, J {/< o «r\ 7|} = o ft" 1 , ij 1 } = [f o Jf} = o. 

Here 1^, ^ff+z+i> -OWl+2; ■■■■> I£n denote the standard basis of the Cartan subal- 
gebra 

TC(f\, / n ), where If = x^i in the case of a hyperbolic summand, If = xf + ^f for an elliptic 
summand, and If 1 = x£i + i — x i+ \^i + T(^i^i + in the complex-hyperbolic case. 

By making a second-order Taylor expansion about I e = I ch = I h = 0, it follows from ( |T3"D 
that Vi = 1, n, there locally exist F y G C°°(£/ ) with {/?, Fy } = {If, F y } = {If 1 , F y } = 
such that: 

(14) /, ° E r >.! ■ 'i ■ E Jr «- i f + E '-w;- 

3=1 j=H+l j=H+L+l 

The non-degeneracy of H(fi, / n ) implies that 

(Fij)(0) G G7(ra; R). 

There is also a parameter-dependent version of this result [El] that is valid in a neighbourhood 
of the orbit, lR n - (vq), and combines the normal form in (|14]) with that in fllCf). More precisely, 
if H.(dpk+i, dp n ) is a Cartan subalgebra, then there exists a neighbourhood, f2, of the orbit, 
]R n ■ (vq), and a canonical map k : Q — » T fc x D with the property that, for all i = 1, n, 

H H+L+l n-k n 

(15) Pi OK- 1 = £iV/j l + £ F ir /f + £ ^./;+ £ i<y • In+l-j- 

3=1 j=H+l j=H+L+l j=n-k+l 

Here I h := (if, i* ), J c/l := Jf +X+1 ) and J e := (J^ +L+2 , I«_ fc ) denote the stan- 

dard generators of the Cartan algebra on E fc , I := ...,Ik) the regular action coordinates 
coming from the normal form in (|10D and D C M. n ~ k a small ball containing G IR n_fc . The 
Fa Poisson-commute with all the action functions. 
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Definition: (i) We say that the rank £ < n orbit, R n • (t> ), is Eliasson non- degenerate 
provided it is compact and the algebra 7i(pe + i, ...,p n ) is Cartan. 

(ii) The rank I < n orbit W 1 ■ (vq) is said to be strongly Eliasson non- degenerate if it is 
Eliasson non-degenerate and in addition, 7r(R™ • (t> )) C M is an embedded ^-dimensional 
submanifold of M. 

2. Quantum integrable systems and Birkhoff normal forms 

Our purpose in this section is to construct a microlocal Birkhoff normal form for a QCI 
(quantum completely integrable) system near a singular orbit. We first recall the definition 
of a QCI system. 

2.1. Quantum integrable systems. Quantum completely integrable Hamiltonians (e.g. 
Laplacians) or those which commute with a maximal family of observables. We now make 
precise the kind of observables we will need to use. 

Given an open U C M™, we say that a(x, £;h) G C°°(U x IR n ) is in the symbol class 
S m > k (U x R"), provided 

We say that a G S™' h (U x IR n ) provided there exists an asymptotic expansion: 

oo 
3=0 

valid for |£| > ^ > with aj(x,^) G S°' k ~ j (U x R n ) on this set. We denote the associated 
H Kohn-Nirenberg quantization by Opn(a), where this operator has Schwartz kernel given 
locally by the formula: 

Op h (a)(x,y) = (2nh)' n [ a(x,£; K) d£. 

By using a partition of unity, one constructs a corresponding class, Op^S 171 ^), of properly- 
supported fi-pseudo differential operators acting on C°°(M). Moreover, this calculus of op- 
erators is independent of the pariticular choice of partition of unity. There exists a natural 
pseudodifferential calculus for such operators with the usual local symbolic composition for- 
mula: Given a G S mi > kl and b G S m2 > k \ the composition Op h (a) o Op h (b) = Op h (c) + 0(h°°) 
in L 2 (M) where locally, 

\a\=0 a ' 

Definition: We say that the operators Qj G Opn(S^' k ); j = 1, n, generate a semiclassical 
quantum completely integrable system if 

[Qi,Qj] = 0; Vl<i,j<n, 

and the respective semiclassical principal symbols qi,...,q n generate a classical integrable 
system with dqi A dqi A • • • A dq n ^ on a dense open set Q C T*M — 0. 
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In this article, we are only concerned with Laplacians on compact Riemannian manifolds. 
Hence, we only consider the homogeneous case where the operators Pi = \fK, P2, P n are 
classical pseudodifferential operators of order one. For fixed b = (b\, b 2 , b n ) G M n , we 
define a ladder of joint eigenvalues of P\ = yA, P2, P n by: 

(16) {(Aifc, \ n k) G Spec(P%, P n ); Vj = 1, .., n, lim = &.,•}, 

where |Afc| := a/ Af fc + ... + A^ fc . We denote the corresponding joint eigenfunctions by (p^. It 
is helpful to rewrite the spectral problem in semiclassical notation. The rescaled operators 
Qj := hPj G Op(S^l x ) clearly generate a semiclassical quantum integrable system in the sense 
of Definition |2~T| . With no loss of generality we may restrict H to the sequence = 1/| Afc|, 
and then the ladder eigenvalue problem (|16|) has the form 

(17) QjtPti — /-^'(^OW; where /i(ft) = 6 + o(l) as ft — >■ 0. 
The semi-classical ladder will be denoted by: 

(18) 

£ b (ft) := {/J,(h) := (jUi(ft), ...,fi n (h)) G Spec(Qi, ...,Q n )] \fij(fi) - bj\ < Cft, j = 1, ...,n}. 

2.2. Model cases. Quantum Birkhoff normal forms are microlocal expressions of a given 
QCI system in terms of certain model system. Model quantum completely integrable systems 
are direct sums of the quadratic Hamiltonians: 

• I h := h(D y y + yD y ) (hyperbolic Hamiltonian), 

• P := h 2 Dy + y 2 , (elliptic Hamiltonian), 

• fch ._ fa [(y 1 D yi +2/2A/2) + y/—l(yiD2 — 2/2A/J] (complex hyperbolic Hamiltonian), 

• I := hDg, (regular Hamiltonian). 

The corresponding model eigenfunctions are: 

• u h (y; A, ft) = \ \ogh\- 1 ' 2 [c + (h)Y(y) \y\- 1 l 2+lX ^l h +c^h)Y(-y) -1/2+^00/*]; |c_(ft)| 2 + 
\c+{h)\ 2 = 1; A(ft) G E. 

• u ch (r,6;t 1 ,t 2 ,h) = |logn|- 1 / 2 r(- 1+i * 1 ^)/ ft e it2 ^ ; t x (h),t 2 (h) G R. 

• u e (y;n,h) = h' 1/A exp(~y 2 /h) $ n (hT 1/2 y); n G N. 

• u reg (9; m, ft) = e imd ; m G Z. 

Here, F(x) denotes the Heaviside function, <$> n {y) the n-th Hermite polynomial and (r,6) 
polar variables in the (2/1,2/2) complex hyperbolic plane. We restrict to those model eigen- 
functions with A(ft), ti(ft), t 2 (ft) = 0(h) as ft — > 0, since these are the ones which localize 
properly. The important part of a model eigenfunctions is its microlocalization to a neigh- 
borhood of x = £ = 0, so we put: 

V>(a;; ft) := 0^( x(aOx(y)x(0 ) • u(y; ft), 

where e > and x G Cg°([— e, e]). In the hyperbolic, complex hyperbolic, elliptic and regular 
cases, we write iph{y] K),ipch{y] h),ip e (y] ft) and ip reg (y;h) respectively. A straightforward 
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computation | T2" | shows that when ti(h),t2(h),nh,mh = 0(h) the model quasimodes are 
L 2 -normalized; that is 

(19) \\Op R (,x(x)x(y)x(0)<y,h)\\v ~ 1 

as h — > 0. Note that, although the model eigenf unctions above are not in general smooth 
functions, the microlocalizations are C°° and supported near the origin. 

2.3. Singular Birkhoff normal form. In this section we introduce Birkhoff normal forms 
for a quantum completely integrable system near a singular orbit. The main result is the 
following quantum analogue (see also |[VN|| ) of the classical Eliasson normal form in (|T5l): 

LEMMA 4. Let c G W 1 be a singular value of the moment map V and ~&L n ■ v be a rank-k 
Eliasson non- degenerate orbit of the joint flow. Then, there exists a microlocally elliptic h- 
Fourier integral operator, F^, and a microlocally invertible nxn matrix of pseudodifferential 
operators, with [Mij,Qk] — 0; k— 1, ...,n and satisfying: 

(20) 

F' 1 (Q 1 - Cl , .., Q n - c n ) F h = n M-(I h - e h {h), I ch - e ch {h), t - e e (h), HD e - e(ft) ) + 0{ff 
Here, I h 3 = hiP, y, ■ y,!),,). J| = tfD% +y], F h = h[(y 3 D yj +y j+ iD y . +1 ) + V^iVjD^ 



y j+1 D y .)) and Ij = KD e .. 

Proof: The proof is essentially the same as in ([VN] Theorem 3.6). The only complication 
here is that since lR n - (vq) is a rank k < n torus and not a point, hDg-^, hDg k must be added 
to the space of model operators. The proof can be reduced to that in [VN] by making Fourier 
series decompositions in the ...,6k) variables (see, for instance |fT2|l Theorem 3). □ 
Consider the microlocal (quasi-)eigenvalue problem in Q : 

(21) Q$ v = n u^ip^j = 1, ...,n. 

We say that Vk{fr) is a quasi-classical eigenvalue if there exists a non-trivial solution of (|2T|). 
The set of quasi-classical eigenvalues around c is thus: 



(22) QT, C {K) := {u(h) : fl2T|) holds, with \u(H) - c| < Ch, j = 1, n}. 

The solution space of (|2l|) can be characterized uniquely (up to a C(ft)-multiple) in terms 
of the model quasimodes ipe^h^ch and ip reg . In the following, we use the abbreviation 
{u e (y; n, h) ■ u h (y; X k (h), H) ■ u ch (y; t 1>k (h),t 2 ,k{K), H) H j=1 e m ^ ) for the expression 

nf =1 ^(%-; h) ® nf+^ chj { yj ; h) ® nf+^f +1 ^(%-; n) ® n n J=L+H+E+l ^ rj {y 3] n) 

in which the (m,n, X,ti,t 2 ) parameters are put in. 

Proposition 5. For any admissible solution i/j u ^ of {\2J\), there exist t 1 (h),t 2 (h), X(h) G R 
and n, m G N and an h-dependent constant c(h) such that 

k 

i/j v = n c(h) F (u e (y; n, H) ■ u h (y; X(H), H) ■ u ch (y; t 1 (h),t 2 (h), h) JJ e im ^ ), 

i=i 

where F h is the microlocally unitary h- Fourier integral operator in Lemma (TJj. 
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Proof. The proposition follows from the uniqueness of microlocal solutions of the model 
eigenfunction equations (see ||CP|| and ||VN|| ) up to C(h) -multiples and Lemma (f|). □ 



3. Blow-up of eigenfunctions attached to singular leaves of the 

Lagrangian fibration 

The purpose of this section is to prove Theorem ([]]) and Lemma ([3]). The key estimate 
is a small scale L 2 mass estimate near a singular leaf A. Roughly speaking, it says that 
the L 2 -mass of the eigenfunctions tp^ with fi(h) G X c (ft) in a tube of radius h s around 7r(A) 
satisfies 

(23) (Op h (xi(x;ti))^, VfM ) > C > 0, 

where C is a constant independent of h and where xi( x \ ^) := Xi(h~ 5x ) with Xi{ x ) a cutoff 
supported near 7r(A). 

In the proof, we will need to use additional pseudo-differential operators belonging to 
a more refined semi-classical calculus, containing cutoffs such as Xi(fr~ Sx )i which involve 
smaller length scales. 

3.1. Eigenfunction mass near non-degenerate, singular orbits. We first give a lower 
bound for the microlocal mass on 'large' length scales of joint eigenfunctions with joint 
eigenvalues fi(h) G £ c (ft) near T sing (c). Let Xreg e C^(T*M) and Xsing £ C™(T*M) be 
cutoff functions supported near r re9 (c) and T sing (c) respectively. We claim that for any 
/ G S(R) with / G Cq°(M) with sufficiently small support and H G (0, Ho], we have 

(24) ^(OpnbCMng)^, ¥V> fi^M^) - c)) > > 0. 

Here, the sum runs over the entire joint spectrum. The estimate (J23|) follows readily from a 
local (singular) Weyl law near a singular rank-fc orbit IR n -t>o. The local Weyl law inturn follows 



from the Eliasson normal form and standard wave trace computations ( ||BPUfl Theorems 
1.1-1.3). When the singular orbit has only hyperbolic summands, the RHS in (|24|) can be 
improved to | logh\ k (see ||BPU|| ). 



Consequently, for each h, there exists at least one joint eigenfunction (p^ of the Q/s with 
joint eigenvalue jj,{h) G T, c (h) satisfying: 

(25) (Oph{Xsmg)<Pi» ¥v)>yr> 0, 

provided we choose the constant C > in (|22|) large enough. By the local Weyl law in 
the regular case, again for C > large enough in fl2~2|) , there also exist joint eigenfunctions 
ipp with n = n(h) G satisfying 

(26) (Oph(Xreg)<P»,<Pn) > ■ 

Here, C7q > denotes possibly different constants from line to line. 
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Definition: Let V sing {h) (resp.V reg (h)) denote the set of L 2 -normalized joint eigenf unctions 
of Qi, ...,Q n with fi{h) G T, c (h) and satisfying fl2"5|) (resp. |26|). 



An easy argument involving invariant measures (see for example |[rZl|| ) shows that the 
connected component T sing (c) must contain a rank- A; < n compact orbit W 1 ■ (v). Again, 
to simplify the writing we can, without loss of generality that there is precisely one such 
compact orbit in T S i ng (c). As the next Lemma shows, the mass much concentrate on at least 
one of these singular orbits A. 

Lemma 6. Let A := W 1 ■ (v ) be an Eliasson non- degenerate orbit, and let Xa( x :0 £ 
C^°(Q; [0, 1]) be a cutoff function supported in an invariant neighbourhood Q of A and iden- 
tically equal to one on a smaller neighbourhood Q of A. Then there exists a constant Co > 
such that for h G (0, ho] and <£> M G V{ft)> 

c 

Proof. With no essential loss of generality, we will assume for notational simplicity that 

V 1 ^) = T sing (c) U T reg (c) 

and that T sing (c) contains a single compact rank-/c orbit, A := R n (i> ). The argument easily 
extends to the case of multiple orbits and components. 



Let Xsing G C™(T*M) be the cutoff functions defined above (cf. |24|- ^5[). Since XiOCsing 



Xa, hence Op h {xk) ° Op h (xsing) = Op h {xk) + 0(h), it follows from (g|) that 

(27) 1/C < (Op h (xa)<P»<Pp) + ((! - Op h (xa)) oOp R (XrtBj)^^/»)- 

Since non-compact orbits do not support probability measures which are invariant under 
the joint flow $ t (see |[rZl|| ), by the compactness of P _1 (c) it follows that A = W 1 ■ (vq), 
the compact rank k < n orbit of vq G "P _1 (c), is the forward limit set for the joint flow on 
7 ?_1 (c). The semiclassical Egorov Theorem then gives: 

((1 - Op h (xA)) ° Op h (xsin g )Vtii <Pn) = (Pp R (ty(l - xa)) ° Op h (x8mg)<Pi» Vv) + 

Moreover, since A = M n • (vq) is an cu-limit set for the joint flow $ t on V~ l (c), there exists 
to G R n such that: 

• supp(&* (1 — Xa)) H supp(xsing) C fi, where the QBNF is valid in Q. 

• supp(%(l - X a)) n R n ■ (v ) = 0. 

In ([T2] Lemma 6), it is shown that microlocal eigenfunctions ip^ given by the QBNF in 
Lemma @ near the orbit A satisfy 

(28) (Opn^Ul - xa)) o Op h (xsin 9 )^, ^) = 0{\ \ogh\- 1 ). 



The lemma then follows from fl27|) , fl28|) . □ 
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3.2. Localization on singular orbits. We claim that the joint eigenfunctions (p^ G V(h) 
must blow up along 7r(A). The first way of quantifying this blowup involves computing the 
asymptotics for the expected values (Op^q)^^,^^) where q G C£°(T*M). 



Lemma 7. Let y? M g V(h). Th 



en: 



(29) (Op n (q)^,^) = \c(h)\ 2 [ qd^ + 0{\ hgh\^ 2 ), 

JM. n -(vo) 

again with \c(h)\ > ^r. 

Proof. Since (p^ solves the equation (|2lj ) exactly (and afortiori microlocally on Q), we may 
express it by Proposition (|5|) in the form: 

(30) y? M = n c(H) F u^, 
for some constant c(H). Here, 

k 

(31) Ufl = (u e (y; n, K) ■ u h (y; X(H), K) ■ u ch (y; h(H),t 2 (H), H) JJ e im ^ ) 

i=i 

Here, by applying the operators on both sides of the QBNF in Lemma @ to the model 
eigendistributions and using the uniqueness result in Lemma ([5]), it follows that for some 
n x n matrix M with M(0) G GL n , 

M(mH,nH,X(H),t 1 (H),t 2 (H)) ■ (mh,nh, X(H),t 1 (h),t 2 (h)) = n{h). 

By the inverse function theorem, the (mh, nH, X(h) , ti(ft) , t 2 {h)) are uniquely determined 
(modulo 0(h°°) ) by the joint eigenvalues fi(h) and moreover, when fi(h) G £(ft) it follows 
that mh,nh, X(h),ti(h),t 2 (h) = 0(h). By Lemma (|), by (|19|) and by (p0|), it follows that 
for ft G (0, Ho], 

(32) ^ < (O^(xa)^,^) = c(h)(F*Op h ( X A)Fu^ Ufl ) < \c(H)\. 

Granted this lower bound on |c(ft)|, the Lemma reduces to estimating matrix elements 
of model eigenfunctions. We now evaluate the matrix elements case by case. The most 
interesting case is where the orbit A is strictly real or complex hyperbolic. Then, given any 
q G C °°(T*M), it follows from @ that 

(33) (Op n (q)<p»,<Pv) = {Opn(q) ° O^(xa)^, Opr(xa)w) + 0(\ log^ 1 ). 

We now use (^Dj) to conjugate to the model setting. The function q goes to go ^ where x 
is the canonical transformation underlying F. The model MJ 1 - action locally reduces to a 
compact torus T fe -action, so we can average the function q o x over the action to obtain a 
smooth invariant function. We then Taylor expand this averaged function in the directions 
(y, if) transverse to the action. We obtain: 

(34) 



(Opn(q) o OpniXn)^^ CWxa) w) = l c (fr) I / qd/j,+ (Op h (r h )u^ tt M ) + {Op h {r ch )u^ u^) + 

i n iv ) 
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where, r h ,r ch G C£°(f2) are the Taylor remainders with r h ,r ch = 0(\y\ + \r)\). A direct 
computation for the model distributions, (see [T2] Lemma 5 and Proposition 3) shows 
that: 

(35) (0 Ph (r h )u„u,) = 0(| log /T 1 /^ (0 Ph (r ch )u„u,) = 0(\ log h\~ l l 2 ). 

The remaining cases are where elliptic (i.e. Hermite factors). Each such factor satisfies 

(36) (Op(r e )w M ,u M ) =0{H) 

so is better than what is claimed. □ 

3.3. Lower bounds on LP norms: Proof of Lemma 0.3. Our objective in this section 
is to refine the argument in Lemma 3.3 to actually produce pointwise lower bounds for 
eigenfunctions attached to Eliasson non-degenerate leaves of the Lagrangian fibration. To 
do this, we study matrix elements on much smaller length scales. That requires the use of 
small scale pseudodifferential operators. We pause to define these objects. 

3.3.1. Small scale semiclassical pseudo-differential calculus. The more refined symbols are 
defined as follows: Given an open set U G R" and < 5 < |, we say that a(x,£;h) G 
Sf{U x R n )if 

(37) \d^a(x,^h)\ <C a ph- 5 ^ + W\ 

Model symbols include cutoffs of the form x{h~ s x, h~ s £,) with x £ C£°(M. 2n ). There is a pseu- 
dodifferential calculus Op?iS™(U x ~R n ) associated with such symbols with the usual symbolic 
composition formula and Calderon-Vaillancourt L 2 -boundedness theorem [Sj]. Composition 
with operators in our original class OpnS m,0 (U x W 1 ) preserves OphS™(U x M n ). 

3.3.2. Outline of proof. Let us now outline the proof. If x^O^j^i ^) e C£°(T*M) is a second 
cutoff supported in a radius h 6 tube around A then clearly 

(38) xl > Xl 

where we view both functions in (|38"D as being defined on T*M and so (|38D is just a pointwise 
estimate. Modulo small errors (see (fEJ)), inequality (P5[) implies the corresponding operator 
bound for the matrix elements: 

(39) ( 0^(xi)w ) > ( Op h (xi)<Pv <Pn )• 

Now, take (f^ G V(h). From the expressions (^) and (0), it follows that the matrix 
elements on the RHS of ( p9[) can be estimated from below by simply computing the masses 
of the model distributions u(y,9;h) = Ilj=i ^"^^(y) • Uh(y) • u c h{y)- Our purpose now 
is to compute these masses on shrinking neighbourhoods of diameter h s centered around a 
singular rank £ < n orbit, A. We show that in such neighbourhoods of diameter h 5 , these 
model distributions have finite mass bounded from below by a positive constant independent 
of h G (0, Ho] provided we choose 5 < 1/2. It follows that, for 5 < 1/2, 

(40) (Opnixl)^^) > ^ >0. 
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Then, combined with (|39|) , the lower bound ( flOf ) implies the L°° lower bound in Lemma 
0.3. We now turn to the details of the proof: 

Proof. By choosing a sufficiently small tubular neighbourhood, M D M inside T*M, we 
can extend the Riemannian metric on M to a Riemannian metric on the entire tube, M. 
Moreover, by possibly rescaling h by an appropriate constant, we can assume that V~\c) C 
M. We will make this assumption without further comment. We will continue to denote by 
A and A e the tubular neighbourhoods defined above. In addition, we define 

(41) A h s = exp{(x, v) G E(5); \v\ < h s }. 

Let xl( x > ^) C^(A h s; [0, 1]) be a cutoff function which is identically equal to 1 in A h s/ 2 . 
Clearly, 

(42) xl(x;h)€S° 5 . 

Now, choose another cutoff function, X2( x 5^ ^)) e C^°(T*M), which is supported in neigh- 
bourhood of radius h s around the torus A = IR n ■ (t> ), with the property that \i — X2 
pointwise. By the Garding inequality, there exists a constant C\ > such that: 

(43) (Opn(xi)^, > (Opnixl)^, ¥V) ~ C^ 26 . 

We now conjugate the right side to the model by the h Fourier integral operator F of 
Lemma (|j). With no loss of generality, we may assume the cutoff function Xii.Ui V> I\ ^) to 
be of product type: 

L+N L+M+N+l n 

xi(y,vJ;fi)=llx(h- 5 y j )x(^%)- U x(h~ 5 pM h ~~ 5 ^) ■ U xi^^n+x-j)- 

3=1 j=L+N+l j= n -l+l 

Here (rj,aj) denote radial variables in the j-th complex hyperbolic summand. Since F is a 
microlocally elliptic ^,-Fourier integral operator associated to a canonical transformation k, 
it follows by Egorov's theorem 

(44) (Opnixi)^,^) = \c(h)\ 2 (Op h (x 6 2 o k)u^) - C^ 25 

where c(h)un(y, 8; K) is the microlocal normal form ( |3~1~D for the eigenfunction (p^. To simplify 
the notation a little, we will write x 5 ( x 'i ^) := x{^~ &x ) below. Now, (x 5 u,u) consists of 
products of four types of terms. The first three are: 



M e = / x \v;K)\x s Ue(v;Wdv, 
j —00 

POO _ 

M h = / x 6 (v,fy\x s Uh(v,?i)\ 2 dr}, 



and finally, 

/oo poo 
/ x 5 (viiV2;h) \x 5 u C h(vi,V2i h)\ 2 d^d^ 
-00 j —00 

To estimate M e , we note that, since cp^ G V(h), and 
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/oo 
e- 2lr > l2/h \$ n (h- 1/2 ri)\ 2 dri + 0(h° 
-oo 



and so for ft G (0, ft ], M e (ft) > 1. 
To estimate Mh, we write: 



(45) M h 



|c + (ft)| 2 +|c + (ft)P 



e -ix x -i/2+i\/h 



X{x/h s £)dx 



log ft y J0 

Since, ip^ G V(H), it follows that |c+(ft)| 2 + |c+(ft)| 2 ^> 1 and the integral in (|45| ) is bounded 
from below by 



(46) 



^0 JO 



0(1 log nr 1 ). 



To estimate this last integral, assume first that £ G [0, ft *]. Then, by an integration by 
parts, 

e~ ix x- 1/2+lX/h dx = o{\h 5 i\ 1/2 ) + o{\n 5 if 2 ), 

'0 

and so, 



(47) 



M h > |logft| 



c x~ 1 / 2+iX / h dx 



odiogftr 1 ). 



'h- s £ 
From (|47|), it follows that: 

(48) M h > |r(l/2 + zA/ft)| 2 (1 - 25) + 0(\ log ft)- 1 ). 

Since, 5 = 1/2 — e, there exists a constant C(e) such that, > C(e) > uniformly for 
ftG (0,ft (e)]. 

Finally, we are left with the integral M c h corresponding to a loxodromic subspace. Since 
\Jk(p)\ < 1 for all k G Z and p G M, it follows that: 



(49) 



M ch > |logft|' 



J k (p)p lt/h dp 



- + odiogftr 1 ; 



Here, Jk{p) denotes the fe-th integral Bessel function of the first kind [AS]. For a > ft 



-d 



(50) 



Mp)p it/h dp\ 



nit/h-p k+i+it/h ^ 



r 



fc+1— it/g\ 



+ 0(|ft 5 ar 1/2 ) =l + 0(|ft <5 a|- 1 



/2> 



and so, 



(51) 



M, 



ch 



/.ft*- 1 , 

> | logft)- 1 / — + 0(| logft)- 1 ) = l-25 + (9(| logft)- 1 ). 

A-* en 



Consequently, given 5 = 1/2 — e it again follows that M c h > C(e) > uniformly for ft G 
(0,ft (e)]. 
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The final step involves estimating (Opft(x a (-0)e lTOe , e im6 ). An integration by parts in the 
I\, I i variables shows that: 

(52) (0 Ph (x\l))e tme , e ime ) = 1 + ©(ft 1 -*). 



As a consequence of the estimates above for My l ,M c y l ,M e and the bound in ([0]) it follows 
that for any e > and 5 = 1/2 — e, there exists a constant C(e) > such that for all 

(53) (Op h (xl)^,^)>C(e)>0. 
Thus, 



(54) \M\U-\ / X d i(x;h)dvol(x)\ >C(e) 



a* 



uniformly for h G (0, h (e)}. Since 



M 



with ft 1 G Spec — yA, the lower bound coming from ( |53D is: 



n — £ 

|^A ? ||L->C(e)A 7 ~' 



The proof of Lemma |3| (i) is complete. Lemma ^ (ii) then follows by applying the Holder 
inequality in the estimate □ 



Remark: From our earlier results in the case of the quantum Euler top [Tl], it seems likely 

n-e n-l 

that in fact Hv^aJIoo > in the elliptic case and ||<^a.,||oo > @^j 4 0-°S^j)~ a f° r an 

appropriate a > in the hyperbolic case. 
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